NACA TN 3298 


DOC. 


yz.H 0 is- ; i 

B'.iSINESS ' _ ___ 

technical >-■ 

NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 3298 

A LOW-DE.NSITY WIND-TUNNEL STUDY OF SHOCK-WAVE STRUCTURE 
AND RELAXATION PHENOMENA IN GASES 
By F. S. Sherman 
University of California 



* 

r 

1 


Washington 
July 1955 


NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 5298 


A LOW- DENSITY WIND-TUNNEL STUDY OF SHOCK- WAVE STRUCTURE 
AND RELAXATION PHENOMENA IN GASES 
By F. S. Sherman 

SUMMARY 


The profiles and thicknesses of normal shock waves of moderate 
strength have been determined experimentally in terms of the variation 
of the equilibrium temperature of an insulated transverse cylinder in 
free-molecule flow. The shock waves were produced in a steady state in 
the jet of a low-density wind tunnel, at initisa Mach numbers of I.72 
and 1.82 in helium and I.78, I.85, 1,90, I.98, 3.70, and 3.91 in air. 

The shock thickness, determined from the maximian slope of the cylinder 

temperature profile, varied from 5 to times the length of the Maxwell 

mean free path in the supersonic stream. A comparison between the exper- 
imental shock profiles and various theoretical predictions leads to the 
tentative conclusions that: (l) The Navier-Stokes equations are adequate 

for the description of the shock transition for initial Mach numbers up 
to 2 , and ( 2 ) the effects of rotational relaxation times in air can be 
accoionted for by the introduction of a "second" or "b ulk " viscosity coef- 
ficient equal to about two— thirds of the ordinary shear viscosity. 


INTRODUCTION 


A precise experimental determination of the structure and thickness 
of the normal shock wave is of fundamental interest in the study of gas 
dynamics because it does much to define the usefulness of the Navier- 
Stokes equations, or of any alternative set of equations, for predicting 
the behavior of a very nonuniform gas. In addition, the observed nature 
of the shock wave sheds considerable light on the magnitude and character 
of so-called "relaxation effects" associated with the finite time required 
to obtain equlpartition of energy among the translational and internal 
motions of a polyatomic molecule. 

The advantages of the normal shock wave as an object of both experi- 
mental and theoretical study are : 

(1) The one -dimensionality of the flow 


2 


NACA TN 3298 


( 2 ) The lack of dependence of the internal shock-wave structure 
upon the nature of boundary conditions at a fluid-solid interface 

( 5 ) The relatively high degree of nonimiformity of the flow^ which 
is sufficient in a shock of moderate strength to cause measurable dif- 
ferences in the profiles predicted by various theories 

(4) The dependence of the degree of nonimiformity, to be defined 
by dimensionless values of the stress and heat flux in the fluid, or by 
the fractional variation of mean flow properties over the length of a 
molecular mean free path, upon a single shock- strength parameter 

It is the last-named property of the shock wave which suggests the 
use of the low-density wind t\mnel for experimental investigation. Since 
the pertinent scale of distance in the shock wave is the mean-free -path 
length, use of the low-density gas stream with its relatively long mean 
free paths will relieve the most troublesome difficulty of shock-wave - 
structure experiments - the accurate resolution of the very small shock 
thicknesses foimd in normal-density wind timnels and shock tubes. 

Whereas the thickness of a shock in air at a Mach nvimber of 2 may be 

1.1 X 10"^ inch at standard temperatiure and pressiire (conditions upstream 
of the wave), the same shock (i.e., same Mach n\miber) will have a thickness 
of 0.0^7 inch in the no. 3 low-density wind tunnel. 

The theoretical problems involved in calculating the profile of the 

steady normal shock wave have been treated in a literatiare^ which orig- 
inated with papers by Rankine, Rayleigh, and Taylor. Subsequently, the 
treatment of the problem by use of the Navier-Stokes equations and the 
assimiptlon of a perfect gas has been broiight to a satisfactory com- 
pletion, so that one may now predict the effects of variation of all the 
significant parameters (Mach number, Prandtl number, and specific -heats 
ratio) and take into accovint the temperature dependence of the viscosity, 
thermal conductivity, and Prandtl number. Even weak relaxation effects 
have been accoxmted for in a recent paper (ref. I 6 ) throiigh the intro- 
duction of a second viscosity coefficient (also called bulk, volume, or 
compression viscosity) which is ass\jmed to have the same temperatvire 
dependence as the ordinary shear viscosity. 

Several authors among those listed have concluded on the basis of 
their calculations that the flow within any but the weakest shocks is 
so nommiform that it cannot be accurately described by the Navier-Stokes 
equations. The equations which are usually suggested as improvements 
are those derived from a kinetic-theory viewpoint, by the methods of 
Enskog, Chapman, Burnett, and Grad (see, e.g., refs. 17 and I 8 ). These 


^See references 1 to I 6 for a representative but not e^daaustive 
bibliography. 
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more complicated equations have been derived for a monatomic gas only 
and have been applied to the calculation of shock- wave profiles in three 
recent papers (refs. 7, 12, and I5) with, results which, were concluded, 
particularly by Grad and ZoUer, to be preferable to the Navier-Stokes 
results for Mach numbers greater than about 1.5» The basis for the pref- 
erence is largely that the predicted shock thicknesses are relatively 
greater according to the more complicated theories. Partly because of 
the plausible natiure of the derivation of the higher order kinetic -theory 
equations, and partly because some of the comparisons of shock thicknesses 
calculated from them and from the Navier-Stokes equations were improperly 
drawn, the Navier-Stokes equations have fallen lately into some disrepute. 
On the other hand, it may be emphasized that this line of thinking has 
been to date almost entirely inspired by theory, since little or no ejiper- 
imental information is available on flows which are sufficiently non- 
uniform to make a measurable difference in the predictions of various 
theories. Even the piirely theoretical trend of thought has encountered 
some rebuttal, notably by Gilbarg and Paolucci (ref. I6). 

Although the usual experimental study of relaxation phenomena in 
polyatomic gases has to do with the absorption and anomalous velocity 

dispersion of forced ultrasonic waves in a gas,^ the effects of lagging 
internal energies are equally evident in the structiare of a normal shock 
wave. If the relaxation times are so small that their effects can be 
represented fluid mechanically by a nonzero value of the b\ilk viscosity 
(refs. 7 and 16), analysis shows that the shock thickness is increased 
by the relaxation effect, quite measurably for even the smallest values 
for bulk viscosity given in the ultrasonic literatxire, although the shape 
of the shock profile is not greatly altered. If, on the other hand, the 
lag of internal energies is so severe that it is appropriate to describe 
the gas within the shock as a reacting mixtirre with sensibly different 
temperatures for the translational and internal motions, then it may be 
expected that the shock profile will assme a shape considerably different 
frcffli that of a shock in a monatomic gas. The profile for the relaxing 
gas woiild be characterized by a zone of comparatively rapid transition 
for the translational motions mixed with and followed by a longer zone of 
adjustment of energy between the translational and internal degrees of 
freedom (refs. 20 to 22). 

The experimental investigation reported herein is not the first to 
be directed toward a determination of the structure and thickness of shock 
waves or to consider the shock wave as a usefiil source of information con- 
cerning relaxation effects. Greene, Cowan, and Homig of Brown University 
have for several years been developing a remarkable program of similar 
Intent, utilizing a measurement of the optical reflectivity of shock fronts 


Reference 19 gives a good summary of this work. 
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passing through a shock tube at atmospheric pressure and above (refs. 23 
to 25). The present studies were vindertaken because the experimental 
method and tools are entirely distinct from those employed by Greene, 
Cowan, and Homig and seem to offer certain advantages in directness and 
detail of measurement and in accuracy and interpretation of resiilts. 

This work was conducted at the University of California under the 
sponsorship and with the financial assistance of the National Advisory 
Committee for Aeronautics. 

s™BOI£ 


Cp 

specific heat at constant pressiore 

Cy 

specific heat at constant volume 

j 

number of excited internal degrees of freedom of a molecule 

L* 

reference length, ft 

Ml 

Mach number upstream of shock wave 

Pr 

Prandtl nimiber 

P 

pressirre, Ib/sq, ft 

To 

wind-tunnel reservoir temperature, 

T* 

reference temperature, ~*'..AlTo 

T 

wire temperature, °R 

tv 

dimensionless, normalized wire temperature 

u 

macroscopic flow velocity, fps 

w 

dimensionless, normalized flow velocity 

X 

distance along flow, ft 

y 

dimensionless distance along flow 
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7 specific-heats ratio, 

5 "area" shock thickness (fig. 10), ft, or units of y 

6m "maximum- slope" shock thickness (fig. lO), ft, or units of y 

6j^j Kronecker delta function 

e dimensionless shock- strength parameter 

K bulk viscosity modulus, Ib-sec/sq. ft; defined by Navier-Stokes 

formulation for viscous stress 


'8u4 8u 


Pij = -P»1J + - f 




.8x-i 8x. 






Maxwell mean free path upstream of shock, ft 

p shear viscosity modixLus, Ib-sec/sq ft 

lx x-component of total molecular velocity, fps 

p density, slugs /cu ft 

Other symbols, used in the appendixes, are defined where they 
appear. 


EXPERIMEINTAL METHOD 


The essence of the experimental method can be described very briefly, 
A normal shock wave is produced in a steady state by introducing an appro- 
priately designed obstacle into the supersonic jet of the low-density 
wind timnel. The shock-wave profile is then recorded in terms of the 
variation in the equilibrium temperature of a small-diameter wire oriented 
parallel to the plane of the shock, as the wire is traversed throiigh the 
shock zone. As a consequence of the relatively long mean free paths char- 
acterizing the low-density wind-tunnel flows, temperature- sensitive wires 
may be obtained with diameters sufficiently small compared with the mean- 
free -path length so that the presence of the wire in the shock zone does 
not disturb the macroscopic flow pattern. Under such conditions of free- 
molecule flow, it is possible to make a theoretical calculation of the 
wire temperature at any point in the gas at which the molecular distri- 
bution function is known and so to relate the wire temperature profile of 
a shock to the velocity or density profile, and so forth. This availa- 
hility of a free-molecvile-flow hot-wire anemometer for low-density wind 
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tunnels was first noted by Tsien (ref. 26), The method is complicated 
experimentally by problems of producing a shock which approximates closely 
the theoretical model, in which the limiting flows upstream and downstream 
of the shock axe perfectly uniform, and by problems of minimizing extra- 
neous thermal exchanges to the wire dvie to radiation and metallic con- 
duction. These are discussed in later sections and in appendix A. 

Theoretical Interpretations of the results stem laxgely from the 
basic work by Stalder and associates (refs, 27 to 29), who investigated 
the free-molecule cylinder in a uniform stream, and from recent work by 
Bell (ref. 50), who considered the effects of nonuiniformities in the flow. 
Appendixes B and C eire devoted to application of their methods to the 
prediction of wire temperature profiles for comparison with the e3q>eri- 
mental results. 


EXPERIMENTAL APPARATUS AND PROCEDURE 
Development of Apparatus 


Shock holders .- The type of shock-producing device which was finally 
developed is shown in figure 1, It consists of a thin-walled circular 
cylinder or cone frustum which is immersed in the uinlform portion of the 
wind-tuinnel jet. The cylinder is alined carefully with its axis parallel 
to the direction of flow and is provided with seme movable device with 
which to vary the open area of its downstream end. When the downstream 
opening is completely plugged, a shock stands detached as in front of an 
impact tube or other blunt body. As shown in figure 2, gradual withdrawal 
of the plugging device decreases the detachment distance and the cuirvature 
of the central portion of the shock wave until a plane shock wave is 
obtained at the entry of the cylinder, (Actually, the shock is swallowed 
in the last picture shown in fig. 2.) The ideal shock holder of this type 
should behave in the following manner; The shock should be held stably 
and exactly across the mouth of the cylinder, so that it will be plane 
and normal to the flow. The shock shouild be followed by a region of utni- 
form subsonic flow of suifficient extent to permit identification of any 
suspected relaxation effects downstream of the shock and then by some sort 
of throat section which will isolate this region from the low-pressure 
region of reexpanded flow behind the shock holder. 

In the development of the shock holder, use was made of the air- 
afterglow flow-visualization technique (ref. 31 ), of a static-pressuure 
tap on the inside of the shock holder near the rear, and of a number of 
very small impact -pres suire tubes. The latter instruments are sensitive 
to the transition of flow properties through the shock because of varying 
viscous and rarefaction effects which are too complicated to be subject 
to theoretical prediction. Since the impact tubes were not small 
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enoiigh to enjoy free -molecule flow, they disturbed the shock wave as 
they passed through, as may be seen in figure 3. The nature of impact- 
pressure traces recorded with a 0.043- inch-diameter tube is shown in 
figure 4, along with a family of ctirves showing the influence of the 
probe position upon the pressure recorded at the shock-holder wall tap. 

As is evident from figiares 2 and 5, different types of choking 
devices were incorporated with the shock holder. The grid shown in 
figure 5 seemed most satisfactory and was used in the final experiments. 

One final item to be noted in the shock-holder design is the slotting 
of the leading edge to permit passage of the measixrlng wire. It appeared 
possible to maJce these slots wide enough to permit good clearance for the 
wires, without disturbing the nature of the shock held across the entry 
plane. 

Limited use was made of another type of shock producer, namely, a 
large transverse cylinder, which produces a detached shock that is locally 
normal in the vicinity of the stagnation streamline but is followed by a 
decelerating flow. 

Free -molecule equillbrlxmi- temperature probe .- Both temperature- 
sensitive resistance wires and butt-welded thermocouples were used to 
fulfill the function of the free-molecule equilibrium-temperature probe 
(fig. 6). The resistance wires have the advantage of availability in 
very s mall diameters and in very strong metals, so that they may enjoy 
free-molecule flow over a wide range of wind-tunnel flow conditions anH 
so tlxat they may be held straight in tension for purposes of alinement 
and location. Since the resistance-wire response is to an integrated 
average temperature over the measuring length, a serious problem of the 
minimization of thermal end losses must be met. This was done in the 
present experiments by introducing separate potential leads to Isolate 
and measiire the resistance of a relatively short central section of the 
current-bearing wire. The potential leads were of the same diameter and 
material (0. 00025-inch tungsten) as the current wire and were very neatly 
joined to the latter in a made -to-order unit furnished by Flow Corporation 
of Cambridge, Massachusetts. The measuring length was about 1/2 inch, 
or approximately the central third of the diameter of the shock-holder 
entry. The potential leads and current wire were oriented in a plane 
parallel to the plane of the shock, and there was no indication that this 
distortion of the probe geometry from that of a simple transverse cyl- 
inder had any effect on the equilibrium temperatijre . 

Thermocouple probes have the advantage of being temperature sensitive 
only at the junction but are not easily available in diameters less than 
0.002 inch because of the mechanical inferiority of most thermocouple 
metals and the diffic\ilty of producing the butt -welded or neatly lap 
welded junction. The range of wind-tunnel operation which will produce 
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free-molecule flow about a wire of this diameter is quite restricted. 
Early models of the thermocouple probe were similar to that shown in 
figure 5« A measuring section 3A inch long, with a junction in the 
center, was held transverse to the flow by its own ends, which extended 
downstream for 1/2 inch before passing into small supporting tubes of 
drawn Pyrex. This arrangement was inferior from a point of view of win 
straightness and of thermal end losses and was abandoned when the feasi. 
bility of slotting the shock holder was discovered. The final thermo- 
couple probe was a simple wire of 0.002- inch iron and constantan, neatl; 
lap welded at the center of a 6- or 7-lnch length (fig. 6(b)). 

Wire traverse .- The final arrangement for support and locomotion 
of the equllibrixmi-temperature probes is shown in fig\rre 7« The ends ol 
the thermocouple wires, or special heavier wires attached to the ends of 
the tungsten resistance thermometer, passed through ceramic insiilators 
mounted on a light spring-steel yoke and were soldered into the measurin 
circuit. The yoke could be spread apart by a differential screw turn- 
buckle, to place the wires in tension, and it was in turn moxmted on a 
carriage which could be run back and forth along the direction of flow 
by a small electric motor. The supposedly constant speed at which this 
traverse moved was 1 inch per 40 minutes, but tests showed that this 
varied slightly over the period of one revolution of the drive motor. 

The position of the wire at any time was given by a dial indicator 
attached rigidly to the upstream wall of the test chamber and bearing 
against the movable base of the wire support yoke. 


Wind-Tunnel Facilities 

The low-density wind tunnel is the no. 3 wind timnel at Berkeley, 
a continuous-flow, nonreturn, open- jet type (ref. 32). The two nozzles 
which were employed for the present tests are designated as nozzle 6 
(1.8 < M < 2.2) and nozzle 8 (3*7 < M< 4.l)(ref. 35)» The gas intake 
to the tunnel is either directly from the room, through activated-alumina 
driers, or from bottled gas. In the latter case, the bottled gas can be 
rendered very dry by running it throijgh a refrigerated trap (ref. 31) at 
about 125 poimds per square inch gage, yielding dew points lower than 
-100° F as measured by an Alnor Type 7OOOL Dew-Pointer. Extra-dry air 
can be obtained by running air from a compressor through this trap. With 
regard to the \iltimate purity of the gases in the wind-ttinnel jet, very 
little of a qiiantitative nature can be said, because of the unknown influ- 
ence of small leaks in the wind-tunnel shell or in the line between the 
wind tunnel and the drier or of outgassing. In general, no special pre- 
cautions were taJcen to insure gas purity diiring the present tests, since 
a few experiments in which nitrogen was kept as pure as possible yielded 
resxilts which were indistinguishable from those obtained with room air. 
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Since the no. 3 wind tinmel is not equipped with special nozzles 
for helium, nozzle 6 was used in this capacity. The flow produced in 
this manner was somewhat less unifom thain the air flow for which the 
nozzle was designed, and the nozzle calibration in helium, is corre- 
spondingly less certain than that in air. 

As mentioned previously, the chief feature of the no. 3 wind tunnel 
which is of interest to the shock-wave study is the production of long 
molecular mean free paths. In figure 8 the mean-free -path length in 
the free stream is plotted, in inches, versus the weight rate of gas flow 
in the two nozzles. Note that the longer mean free paths occur in the 
lower Mach number nozzle, which fact made it easier to obtain shock- wave 
traces in nozzle 6 and accounts somewhat for the scanty data obtained 
for the higher Mach number flows. Another factor in the wind-tunnel per- 
formance which limited the scope of tests in the higher Mach number nozzle 
was a slight fluctuation in the test-section pressure which became serious 
at flow rates greater than 10 poimds per hour, making it impossible to 
produce a perfectly steady shock. Something of this nature occurred also 
in nozzle 6 but seemed to be the fault of the shock holder, whereas even 
a detached shock off a transverse cylinder was -unsteady at high flow rates 
in nozzle 8 . 


Technique and Schedule of Final Experiments 

After a fairly extended period of preliminary and developmental 
experiments, a series of final tests was performed. The testing procedure 
•was very simple once a shock holder and thermocouple or resistance-wire 
probe had been installed and alined. This alinement was easily achieved 
when a thermocouple was used, since the shock holder was mounted in a 
maimer allowing delicate adjustment of its height and inclination and 
since there were onky two slots in the shock holder for the wire to clear. 
The resistance-wire-thermometer alinement, requiring clearances in four 
slots, was somewhat more tedious and was accomplished by adjusting both 
the shock-holder orientation and the positions of the ends of the fine 
tungsten wires. The latter adjustment was achieved by bending the hea-vy 
wires which were attached to each tail of the tungsten wires. 

When the wind tunnel was ready to be operated, the free-molecule cyl- 
inder was first placed in the entry plane of the shock holder to provide 
an indication of proper location of the shock -wave. With flow established 
at a Mach number and pressure level known to be appropriate to the par- 
ticular shock holder in use, the shock-holder choking grid would be set to 
give first a detached shock and then a swallowed one. The wire-temperature 
indications for these two conditions corresponded approximately to the 
do-wnstream and upstream limits of the shock-wave profile, and, when the 
choking grid was further manipulated to make the temperature of the wire 
(still in the entry plane) nearly equal to the average of these end 
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readings, the shock was considered to he properly placed. One factor 
due to the influence which the shock position had on the balance between 
jet static and test -chamber pressures complicated this procedure. Since 
the shock holder captured a sizable fraction of the uniform portion of 
the jet, there was a considerable interaction between the shock position 
and the aforementioned pressure balance, making it necessary to adjust 
the two more or less slmtataneously. 

Once the shock had been positioned, the wire temperatiire was observed 
over a period of time, in order to make certain that the shock was sta- 
tionary. Then the measuring wire was run downstream until a definite 
indication of the uniform flow downstream of the shock was reached. From 
this point, the shock profile was traced out in the upstream direction by 
two different techniques, one appropriate to the resistance thermometer 
and the other, to the thermocouple. 

Shock traces with a resistance thermometer .- The resistance wire was 
fed a constant ctirrent of 0.2 milliampere, obtained from a 5-volt dry 
cell, a fixed 10,000-ohm resistor, and a 0- to 12,500-ohm variable resis- 
tor in series. The potential drop across the temperature-meas\rring sec- 
tion of the wire was measured with a Rubicon hand-balancing potentiometer 
and sensitive galvanometer (Leeds and Northrup, 0.44 microvolt per 
millimeter), while the wire current was monitored by a Weston Model 31 
milliammeter. The potentiometer readings were converted to temperatures 
by use of a linear formula with a temperature coefficient of resistance 
which was determined by the wire manufacturer from calibrations performed 
on samples from the same spool of wire. According to this calibration 
the least count of the Rubicon corresponded to about 0 . 1 ° F when the wire 
current was 0.2 milliampere. 

While the wire temperature was being measured in the above manner, 
the stagnation temperature of the flow was measured by an iron-constantan 
thermocouple located in the upstream reservoir chamber. Before each run, 
an initial reading of wire resistance versus thermocouple temperatxare was 
taken at no-flow conditions, with a siifficiently high gas pressure to 
make the wire temperature correspond to the true gas temperature. The 
scatter of these no-flow readings from day to day gave some idea of the 
probable error of determination of the absolute temperature of the wire 
due to calibration drifts, instrument errors, and so forth. 

During the tracing of a shock profile, the potential drop across the 
resistance wire and the wire position were recorded at a sufficient number 
of points, the spacing of which was determined by the mean-free -path 
length, to give a satisfactory definition of the shock profile and thick- 
ness, as shown in figures 9(a) to 9 (h). The time reqiiired for this was 
about an ho\ir per shock wave. 
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Traces vith thermocouple vire .- When thermocouple probes were used, 
the difference in thermal electromotive forces generated by the probe 
jimction and by an identical junction placed in the stagnation chamber 
was automatically recorded by a Brown Electronik X-Y function plotter. 

The recorder chart was driven at constant speed by a small motor wired 
in parallel with the motor on the wire traverse. Distance on the chart 
was converted to distance traveled by the wire by assuming wire and chart 
speeds to be constant at values measured with a stop watch. Each shock- 
wave trace was recorded in both the ij^jstream and downstream directions, 
and when the recorder sensitivity was adjusted to the maximum limit for 
stable pen motion the hysteresis between forward and rearward running 
traces amounted to a few thousandths of an inch of wire travel. The 
stagnation temperature was recorded separately with the same thermocouple 
as was used with the resistance wire. 

Schedule of final experiments . - The final data which are reported 
herein were taken at initial (or upstream) Mach numbers of 1 .J 2 and 1.82 
in heliiam and I.78, I.85, I.90, I.98, 3-70, and 5.91 in air. 


The dimensionless variables introduced for analysis of the data were 
suggested for the most part by Grad's analysis (ref. 15 )• The distance 
in the direction of the flow is converted to a variable y, defined by 


This variable is discussed in appendix B. In the formula for y, 

X is the measixred distance, L* is a reference length having the same 
units as x, and e is a parameter indicating the strength of the shock; 
L* and e are computed from more familiar variables by the equations 


REDUCTION OF DATA 





( j^ + 8j + 15 )(Mi^ - 1) 
(j + 3 ) + (J + 5 )Mi^ 


e 
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The gas properties which have been assumed for the reduction of the data 
are as follows: For helim: 

j = 0 , Pr = 2/3, K = 0 

and for air: 

j = 2 , Pr = 5 A, t = 2 p /3 

The temperature dependence of |i was taJcen from references 54 and 35. 
Table I gives a listing of the calculated values of y/x for each flow 
setting at which a final shock profile was traced, along with other var- 
iables needed to describe the flow. 

If one prefers the upstream mean free path Ap as a reference 
length, he may convert the abscissa of the shock profiles y to x/Ap 
Piy dividing y by y/x and then by Ap. Both of the last-named (quan- 
tities are listed for each flow setting in table I. 

The free-molecule cylinder temperat\rre is determined from the cali- 
bration of the thermocouple or resistance wire and is made dimensionless 
by forming its ratio to the measxared stagnation temperatiire . This ratio 
is then normalized to run from 1 to - 1 , in the manner applied to the 
theoretical profiles in appendix B, by use of the formula 



in which subscripts refer to conditions at the upstream (l) and down- 
stream (2) limits of the shock zone. This variable artificially cancels 
the effects of calibration differences and radiation and end losses, as 
well as the theoretical variation in temperature spread across the shock, 
simplifying comparisons of experiment with theory or of one experiment 
with another. Introduction of the variable is in keeping with the aim 
of the present investigation, since the normalization does not influence 
those characteristics of the shock profile (i.e., shape and thickness) 
which are of principal interest herein. The normalized temperature can 
be formed only for the traces of normal shocks produced in the hollow 
shock holder, because the detached shock traces have no point of zero 
slope to serve as a reference level on the downstream side.. 

The normal- shock traces are presented as graphs of tv versus y, 
based upon the step-by-step data of the resistance-wire measurements and 
upon a curve faired smoothly through the middle of the hysteresis loop of 
the thermocouple traces. In drawing the latter curve, an attempt was made 
to remove the effects of the small periodic speed variation in the wire- 
traverse motor. The origin for the variable y is taken arbitrarily at 
about the point at which tv eq.uals zero. 
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The shock thicknesses are determined graphically from the curves 
of t^f versus y by both the ma.xlm\jm~slope method (Prandtl's definition) 
and the area method recently suggested by Grad. Figure 10 compares the 
two definitions, which yield the same numerical values only if the shock 
profile is a straight line. For any other profile the area thickness is 
greater than the maximm- slope thickness; for example, for the zero- 
strength profile (see appendix B) the area thickness is 2.77 units 
of y, while the maximum-slope thickness is 2. In what follows, the area 
thickness will be designated by 6 and the maximum- slope thickness, 
ty 6m. 


EXPERIMEKTAL RESULTS 

Free-Molecule Probe Response In a Uniform Stream 


Both theory and previous experiment (refs. 27 to 29) have shown 
that a cylinder which Is a perfect heat conductor Internally but Is 
totally Insiilated from radiation and end losses. If placed with Its axis 
perpendicular to a uniform stream of gas In which the mean-free -path 
length Is several times the cylinder diameter, is heated by the stream 
to a temperatiire which Is a function only of the local Mach number and 
static temperature and of the number of excited degrees of freedom of the 
gas molecules. For a stream with constant stagnation temperature, the 
cylinder temperature increases fairly rapidly with Mach number up to 
about M = 2, after which the response is q\aite insensitive to further 
increase in M. The first Important results of the present experiments 
demonstrate the extent to which this behavior was found in the wires 
which were used as shock-wave probes. These wires were different from 
those with which the free -molecule theory was tested by Stalder, Goodwin, 
and Creager (ref. 29) in that no attempt was made to combat thermal end 
losses by artificial end temperatttre control. Aside from the special 
arrangement of the resistance-wire thermometer, with extra potential leads 
to isolate a central portion of the ciarrent wire, the present wires 
depended only on a high length-diameter ratio to reduce end-loss effects. 
The results obtained axe seen in figure 11, in which the ratio of measured 
wire temperature to measured stagnation temperatTore is ccBnpared with the 
theoretical ratio for a perfectly insialated cylinder. This figure shows 
n.l 1 experimental points for which the Knudsen number was greater than 
about 5 lying as close to the theoretical curves as can be expected in 
view of probable errors in measurement. (Exceptions to this statement, 
for which no explanation has been found, exist in the case of certain 
data for helim at supersonic Mach numbers.) As the Knudsen number drops 
below 5^ the points tend to drop below the theoretical curve because of 
a lack of fully developed free-molecule flow. 
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As a consequence of tlie results shown in figure 11 it is thought 
that the wires used in this experiment behaved as insulated cylinders 
in fully developed free-molecule flow, for a1 1 wind-tunnel conditions 
encoimtered, as long as the mean free path in the vicinity of the wire 
was at least 5 wire diameters long. 


Operation of Normal Shock Holders 

Experimental difficulties which were encoimtered in the attempt to 
realize the ideal shock-holder behavior described earlier appeared to 
be of two kinds, both probably traceable to the boundary- layer growth 
inside the shock holder. First, if the internal -boimdary -layer b-uild-up 
was excessive, the shock holder became self -choking and it was Impossible 
to draw the detached shock back into the entry plane. Some allowance for 
this boundary-layer growth was made by diverging the internal contour, 
but no veiy satisfactory theory exists to guide the correction, which at 
best woiild be applicable over only a very small range of flow conditions. 
Second, there existed certain conditions, possibly corresponding to an 
overcorrection for boundary layer, for which the shock woiold not remain 
stationary in the shock-holder entry but performed small random oscil- 
lations about a mean position. 

Fortunately, almost every shock holder which was tried worked prop- 
erly over some range of wind-tunnel conditions, so that by use of the 
three models shown in figure 1 it was possible to cover an adeq-uate over- 
all range of flows. However, even when a shock holder was operating cor- 
rectly, with the shock steady in the entry plane, the internal boundary 
layer exerted a slight influence on the downstream tail of the shock-wave 
profile, as will be disctissed more fully later in this report. 


Normal-Shock-Wave Profiles 

The normal-shock traces in raw-data form showed wire-tenperature 
changes of about 60° F across the shock and shock thicknesses varying 
from a few hundredths to a few tenths of an inch. The final traces are 
presented in dimensionless, normalized form, as plots of t-yj- versus y 
(see "Reduction of Data"), in figures 9(a) to 9(h). Each figure shows 
all traces recorded at a given flow setting, including both thermocouple 
and resistance-wire traces and occasional repetitions of the same trace 
which were taken on different days or with different shock holders. The 
latter repetitions were frequent in the higher Mach number nozzle, because 
of the large amount of scatter which was encountered. For point -by-point 
traces with the resistance wire, the points of measvirement are shown; for 
the continuous thermocouple traces recorded with the X-Y function plotter, 
the smooth curve is shown. (An exception to the last statement is seen 
in fig. 15 (a), where points have been picked off the thermocouple trace 
to avoid confusion with the theoretical curve . ) 
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Shock-Wave Thicknesses 

Shock thicknesses were calciilated only for those traces which were 
obtained with the resistance wire, since the higher Khudsen n\jmbers 
obtainable with this wire gave greater assurance that the wire did not 
perturb the natural condition of the shock front. These results are 
shown in figure 12 , in which 6 is plotted versus the shock- strength 
parameter e; in figure I5, in which ^i^/pubm is plotted versus the 
inltiail Mach number Mi; and in figure 14 , in which 5 m/-'^l is plotted 
versus M]_, In all graphs the results of repeated measurements are rep- 
resented by a point giving the average value and a line showing the 
spread of values. Various theoretical results from the literature, and 
experimental results by Greene, Cowan, and Homig, are also shown in 
figure 15. In examining figure 15, one should note that some of the 
shock thicknesses shown are taken from the density profile; some, from 
the velocity profile; and some, from the free-molecule wire temperature 
profile. The effect of these Inconsistencies is discussed somewhat in 
a later section. 


Evidence of Wire-Size Effect on Shock-Wave Profile 

The range of Khudsen numbers (h/d) encountered in the present 
experiments was sufficient to indicate in some instances the nature of 
wire-size effects when the conditions for fully developed free-molecule 
flow were not met. Such effects may be seen by contrasting figure 9 (a)^ 
which shows identical traces obtained with a relatively large thermo- 
couple wire (A*/d = lo) and with a much smaller resistance wire 
^A*/d = 80^ at a flow condition for which both wires enjoy free-molecule 
flow, with figure 9(g), wherein the trace from the larger wire corre- 
sponds to A^^d = 1.5 aJid that from the smaller wire corresponds to 
A*y d = 11. Evidently the larger wire in the latter case has produced 
an artificial broadening of the shock front. Such comparisons between 
traces of identical shocks, recorded with wires of different sizes, in 
addition to the study of deviations between experimental and theoretical 
wire temperatures in the uniform stream, served to establish the Khudsen 
number criterion for discarding any shock traces which could contain a 
wire- size effect. 

EXPERIMENTAL ERROR 
Random Error in Measurements 

The three quantities measured were two ten5>eratures and a distance. 
In addition, use was made of the previously determined calibrations of 
the wind-tunnel jets. 
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The accijracy of determining the absolute level of ten5>erature, with 
both thermocouples and resistance wires, was quite poor, with errors 
amounting to perhaps ±5° F. These errors remained approximately constant 
ng a run and were mostly canceled by the method of data reduction. 

The sensitivity of the temperattire-measviring equipment, on the other hand, 
was quite good, being about ±0.1° F for the resistance wire according to 
the least coimts of the potentiometer and galvanometer, which corresponds 
to about 0.16 percent of the total change in wire temperature across a 
shock. Thp stagnation-temperature measurement also probably contained 
ei*rors of a few degrees, but again the data— reduction method involves only 
changes in Tq during the tracing of a shockj and these changes could be 
measured with a precision of about ±0.5 F. 

Distance through the shock was measured with an effective least count 
of about 0.0002 inch at the higher Mach numbers (O. 0001- inch- reading dial 
indicator) and of about 0.001 inch throughout the remainder of the tests. 
The smallest change in x which could be reliably measured amounted to 
2 percent or less of the maximum-slope shock thickness in all cases. 

Uncertainty in the calibration of the wind-tunnel jets entered into 
the conversion of x into y and into the specification of the initial 
Mach number. Previous experience with nozzles 6 and 8 indicates that 
there exists a probable error of ±2 or ±3 percent and of less than 1 per- 
cent in their respective Mach number determinations. Corresponding to 
this Mach number uncertainty and to the probable error of measurement in 
stagnation pressirre and temperature, there are probable errors in L* 
onti A-ji^ of about ±5 percent in nozzle 6 and of about ±1 percent in 

nozzle 8. The probable error in the conversion factor y/x is consid- 
erably less, being less than 1 percent throughout, because of the fact 
that y/x- is practically independent of M]_ within the range of nozzle 6. 

The greatest source of uncertainty in the experimentally determined 
shock thicknesses entered through the difficulty of drawing the maximum- 
slope line. This error appeared to be about ±5 percent up to about Mach 
number 2 and ±10 percent at the two points around Mach n-umber 3.8. 


Systematic Experimental Error 

Any new e3q>erlmental method must be scrutinized carefully for sources 
of systematic error, and the number of possible sources of such error is 
great in the present experiment. It may be that the following discussion 
will fail to riamp them all, but everything which has occurred to the author 
individually or has been suiggested by others will be given some mention, 
even if some items are concluded to be unimportant within the small range 
of tests performed so far. This discussion will be presented as an attempt 
to answer a nvtmber of critical questions concerning the experiment. 
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(a) Was the shock vave really plane and normal to the incident flov? 
Since "both the thermocouple and the resistance wires recorded the result 
of an average energy balance over an appreciable length of wire, the 
experimental tv traces would certainly be distorted from the true shock 
profile if the shock were curved in such a manner that the wire did not 
lie parallel to the shock front. By reason of this same averaging prop- 
erty, the wire itself could not be used effectively to check the plane- 
ness of the shock (for instance, by varying the position of the thermo- 
couple junction along the length of the wire). The principal assurances 
that the shock was truly plane in the present tests were as follows: 

(1) The shock stood directly in the entry plane of the hollow shock 
holder, according to measurements of the wire temperature. The shock is 
known to be convex upstream when slightly detached, and it appears rea- 
sonable that it woiold be convex downstream when very slightly swallowed; 
the inference is that it would be plane when in the entry plane, (2) Com- 
parison with the profiles of detached shocks from the transverse cylinder, 
for cases in which the detachment distance was fairly large compared with 
the shock thickness, showed good agreement, although what shock curvatures 
there could have been in the two cases were certainly very different in 
nature. Neither of these assinrances amounts to a direct observation, and 
the planeness of the shock may still be open to question; in any case it 
is evident that the error introduced into the shock-wave profiles would 
be an artificial thickening. 

(b) Was the shock really steady, or could it have been drifting or 
vibrating? High-speed schlieren photographs show that shocks in 
atmospheric -pres sure wind tunnels frequently perform rapid oscillations 
about a mean position, so that shocks appear considerably sharper in 
high-speed photographs than in low-speed ones. At the other extreme, 
shocks have been observed to drift slowly in one direction in the low- 
density wind tunnel, in response to a transient change in the very thick 
boundary layer during the time a model requires to reach its equllibriimi 
temperature in the supersonic stream. Either one of these effects, if 
present during the tracing of the shock profile, wovild raise serious 
questions concerning the validity of the data. Errors due to a slow drift 
could be avoided by waiting for the transient to pass, and considerable 
care was taken to do this. Vibrations, if fairly slow, would show up on 
the wire temperature, particiilarly when the wire was near the center of 
the shock. Such slow vibrations could be caused by mechanical vibration 
of the shock holder or of the choking device, or by random fluctuations 

in the test-chamber pressure. In either case they were easy to identify 
and avoid. High-frequency vibrations, on the other hand, coxild not be 
detected by the rather sluggish measuring circuits, and the. only argu- 
ments against their possible existence are indirect. No one has as yet 
sxiggested any reasonable source of excitation of such high-frequency 
oscillations (the low-density supersonic jet is almost certainly not 
turbulent), and in any case the error Introduced would again be a broad- 
ening of the shock front, a resTJit which appears somewhat unlikely in 
view of the rather small shock thicknesses which were observed. 
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(c) Was the measured shock profile unduly influenced hy the houndary- 
layer growth within the shock holder or hy the geometiy of the shock 
holder itself? During the coijrse of preliminary and developmental experi- 
mentation it became evident that the shock was almost never followed hy a 
long zone of perfectly uniform subsonic flow but rather hy locally accel- 
erating or deceleratiixg flows, the exact nature of which depended upon 
the shock-holder geometry and the Reynolds number. Furthermore, at the 
lowest Reynolds numbers the effect of these nonuniform downstream boundary 
conditions appeared to penetrate upstream to a considerable distance 
through the shock wave, as evidenced by comparisons of normal shocks and 
detached shocks. However, final data were gathered only under conditions 
of mirrlmurn downstream nonuniformity; and comparisons between normal- shock 
traces recorded with shock holders of different sizes and shapes, or with 
the same shock holder at different Reynolds mmibers, suggest that this 
effect is not very large. Nevertheless, this is one of the most real 
disadvantages of the present experiments and one which woiiLd merit con- 
siderable attention in an attempt to refine or extend this work. 

(d) Coiold the tv profile have been significantly distorted by 
x-dependent values of radiation and end conduction losses or of Khudsen 
numbers? Since the wire temperature, effective end temperature, radiant 
surroijnding temperature, and convective-heat-transfer coefficient all 
changed as the wire was moved through the shock, partlciilarly as the wire 
first passed the shock-holder entry plane, changing amounts of extraneous 
heat losses coilLd conceivably have pertvirbed the shape of the shock pro- 
file. However, some direct calciilatlons (see appendix A) and a nmber 

of indirect observations indicate that this effect was unimportant. The 
indirect evidence which is considered pertinent is as follows: (l) Gtood 

correlation between shock profiles at various gas pressures, for which 
the calculated effects of extraneous heat losses were significantly dif- 
ferent and for which the distances of wire travel were also distinctly 
different; (2) lack of any discontinuity of profile slope or curvatvire, 
to be associated with the wire's passing the entry plane; and (5) good 
agreement between thermocouple and re si stance -wire traces, end losses 
from the two instnmients being different in calculated magnitude. Effects 
of varying Khudsen number may also be disregarded, since final data were 
retained only for traces in which the Khudsen number was everywhere above 
the free-molectile flow limit. 

(e) Was the shock profile unduly influenced by local nonuniformity 
in the free-stream flow? This question has particular significance with 
respect to the use of an air nozzle to produce helium flows, which were 
not so uniform as might be desired. However, tests made with the shock 
holder in different axieil positions (and hence in different situations 
regarding local nonuniformities of flow) produced essentially identical 
results, so that this possible sovirce of error may apparently be dismissed. 
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DISCUSSION 

Significance of ty Profiles 


In assessing the value of the present experimental method, one is 
bound to question the relative usefulness of knowing the t^^ profile of 
a shock wave as contrasted with knowing the usual profiles such as those 
of velocity or temperature. This is obviously a difficult point to dis- 
cuss conclxosively, and the suggestions which follow cannot be expected 
to satisfy all critics; nevertheless, the point is too important to be 
passed by without mention. 

There is at first the question of tv being, for a cylinder of 
finite diameter, as good a point fiinction as (say) the velocity. A theo- 
retical argument here is difficult, since one is not frequently critical 
in defining how satisfactory a point function the velocity itself is; 
but a practical, answer may be gained experimentally by observing the 
change in tv which occurs when a given wire is displaced by 1 diameter 
across the steepest gradients in the shock. If this change is signif- 
icant, then the wire is presumably too large to yield a satisfactory 
point measurement. (The largest such change was, at M-, = 3.9, about 
1 

2— percent of the total change across the shock, a value which may already 

be excessive. The next largest was about 1 percent, at = 3.7, »nd 
all others were considerably less than 1 percent.) 

Supposing that tv is a satisfactory point function, one then asks 
how much information about the complete structure of the shock wave is 
gained by knowledge of the tv profile. At the outset, it is evident 
from the calculations outlined in appendix B that there does not exist, 
except in the limit of vanishing shock strength, any inversion relation 
which will yield the velocity or temperatinre profile directly from the 
tv profile. From the viewpoint of continuimi theory there appear to be 
two reasons for this fact, the first being the general lack of locally 
adiabatic flow within the shock and the second being that the wire tem- 
perature is influenced by local values of the stress and heat flux in 
addition to the velocity and temperature. On the other hand, the cal- 
culations have shown that the tv profile does not differ greatly from 
the w profile for air or helium at a Mach number up to 2 , regardless 
of the particular theory used for the calciolation. Some distinguishable 
differences in profile shapes for tv and w arise from the thirteen- 
moment theory at a Mach number of I.6I and from the Navier-Stokes theory 
for Mjl = 3*7# t)Ut under any circmstances it appeeirs clear that the wire 

temperatvire is more akin to the gas velocity t han to any other variable. 
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A more general viewpoint which is subscribed to herein is that the 
state of the gas at a point within the shock is determined when the molec- 
ular distribution function is known at that point and that, of all the 
observable mean properties (such as velocity, density, gas temperat\ire, 
or free-molecule wire temperatiare ) which can be calculated from tte dis- 
tribution function, no single one plays a role essentially more signif- 
icant than that of any other. The function of an experiment in which 
only one of these observable quantities is actually measured is to infer 
the correctness of some theoretically postulated distribution function, 
in the hope that a distribution function which predicts the correct valxies 
for the measured observable quantity will also predict correctly those 
observable quantities which coilLd not be measin-ed. Th-us, if it is said 
in this report, for example, that the Navier-Stokes equations appear sat- 
isfactory to predict the shock transition up to a Mach number of 2 , what 
is meant is that it has been possible to find a distribution function 
which yields the Navier-Stokes equations when substituted into the Maxwell 
transport equations for mass, moment\mi, and energy and which, when used 
in the calculation of the free-molecule cylinder temperat^n•e, yields 
results in good agreement with experiment up to a Mach number of 2 , The 
inference then made is that the velocity and temperature profiles which 
are calculated en route to the successful prediction of the tv profile 
(see appendix B) axe also correct. 

The danger of this line of thinking lies in the obvious lack of 
unique correlation between a correct distribution function and a correct 
derived value for any single observable quantity, particularly when the 
range of experimental data is sharply limited, as it is in the present 
case. In other words, it may be possible to find two or more distinctly 
different distribution functions from which one could calciilate practi- 
cally identical tv profiles for a shock in a given range of Mach number, 
whereas these same two distribution functions might lead to vastly dif- 
ferent (say) temperature profiles for this shock. Fortixnately, no such 
examples have appeared as a result of the theoretical calculations done 
so fax. 

The negative corollary of the above viewpoint may be a little more 
readily acceptable,* this is, a distribution function which leads to a 
very poor prediction of the observed tv profile may probably be dis- 
carded as wrong. Even this conclusion must be tempered by one's estimate 
of the adequacy of the technique by which the parameters of the distri- 
bution function are determined (by use of the Maxwell transport eq^ua- 
tions, etc.). 

The most positive suggestion which seems appropriate for this report, 
in view of the presently immature condition of both the theories and 
experiments, is that the t^ profiles may be used to infer the general 
superiority of one theoretical approach over one or more alternatives and 
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that the free-moleciile wire- temperature measurement is about as useful 
as a velocity or density meas\arement for the purpose of this inference, 
at least for Mach numbers up to about 2. 


Comparison of Experiment With Theory 

As noted in appendix B, theoretical shock profiles for comparison 
with experiment have been calculated for = 1.82 in heliimi, = I.89 

in air, and M]_ = 5.7O in air. The calciilations were done by the Navier- 

Stokes and Mott -Smith methods, as indicated in table II. Figures 15 (a) 
to 15(c) show graphical comparisons of the theoretical ty^ profiles with 
the experimental shock traces at corresponding Mach nmbers. For the 
helium shock, and for air at Mj = I.89, a veiy satisfactory fit to the 

data is afforded by the Navier-Stokes theory, with zero bulk viscosity 
in the monatomic gas and with bulk viscosity eqiial to about two-thirds of 
shear viscosity in air. The shape of the Mott-Smith profiles is satis- 
factory, but the shock thickness is evidently too large, at least for 

the transport equation (see appendix B). At M-|_ = 5*70, none of 

the theories yields satisfactory agreement with the data; the Navier- 
Stokes profile has the wrong shape and the Mott-Smith profile is still 
too thick. 

The thirteen-moment eqiiations have not been solved for a normal shock 
wave for Mach numbers greater than I.65, a limit which is unfortimately 
lower than the lowest Mach number obtained in the experiments. In order 
to obtain some idea as to its correctness, this theory and the Navier- 
Stokes theory are compared for helium at Mj_ = I.61 in figure I6. At 

this Mach number, the shape and thickness of the thirteen-moment profile 
are distinctly different from what would be expected from an esctenslon 
of the experimental data toward lesser shock strength, assuming that good 
agreement between experiment and the Navier-Stokes theory continues in 
this region. 


Relaxation Phenomena 

Bulk viscosity .- The ciorves displayed in figure 15(b) strongly sug- 
gest that air shoiold be characterized by a bulk viscosity coefficient as 
well as a shear viscosity. Although the concept of a nonzero bulk vis- 
cosity coefficient is frequently overlooked in the study of fliiid dynamics, 

3 

it has been the subject of many papers, both theoretical and experimental.^ 


^Thls literature has been very well summarized by Truesdell (ref. 58 ; 
see sec. 6 IA, pp. 228-251). 
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The conclusion to be drawn from this literature is that bulk viscosities 
certainly do exist in any but monatomic gases but that this fact is of 
little practical importance in all but a few flow situations. The excep- 
tional circumstances are those that have been noted before: High-frequency 

sovind absorption and dispersion, shock waves, or any flow in which the 
normal, rather than tangential, viscous stresses are predominant. Exper- 
imental values of bulk: viscosities, as obtained from acoustical interfer- 
ometry and related experiments, have been very recently collected and 
tabulated by Truesdell (ref. 19 ). The value given for air is k = 0 . 57 ^, 
in good agreement with the value obtained by the present study of shock- 
wave structure. 

According to several writers (e.g., Wang Chang and Uhlenbeck, 
ref. 57), the ratio it/n can be related to the relaxation time associ- 
ated with the internal degrees of freedom when these degrees of freedom 
exchange energy "easily" with the translational degrees of freedom of the 
molecule (i.e., when relaxation times are small). From the relationships 
given in reference 37, one calctlates that a bulk viscosity k = 2/5 d 
implies the average requirement of 5*3 molecular encovinters for the reduc- 
tion of a sudden disturbance in the rotational molecular degrees of free- 
dom, in air, to l/e of its initial magnitude. 

Additional relaxation effects .- Althovigh the present findings con- 
cerning the importance of bulk viscosity in the shock wave are corroborated 
by the results of Greene, Cowan, and Hornig's optical- reflectivity exper- 
iments, the latter experiments seem also to have suggested some additional 
or residual relaxation effects downstream of the main transition zone. 

This has led to the feeling that the shock transition in certain diatomic 
and polyatomic gases may occur as a two-stage process, in which 90 percent 
or so of the total density change occinrs relatively quickly as described 
by the Navier-Stokes equations with bulk viscosity, after which some final, 
slower adjustments are completed to bring about the final equilibrium 
state. Although the appearance of the wire-temperat\ire traces would seem 
to deny the existence of such a two-stage shock transition in air or 
nitrogen, certain featiares of the present experiment suggest that it may 
be incapable of showing up the small adjustments of the second, slower 
stage, even if a two-stage process does exist. An interesting point is 
brought out by reference to appendix C and to figure 17, which suggest 
the possibility that the wire temperatiore in the downstream tail of the 
shock may be almost independent of the state of relaxation of the gas, 
particularly in the Mach nmber range in which most of the present data 
were recorded. (See particularly curves a and c in fig. I7.) While the 
results of appendix C are certainly Inconclusive, because of the roiaghness 
of the relaxation model and the lack of knowledge of values for the var- 
ious accommodation coefficients, there are two other features of the pres- 
ent work, namely, the limited experimental accuracy and the distortion of 
the downstream tail of the shock by boundary- layer growth in the shock 
holder, which would definitely hinder detection of a second stage in the 
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shock transition. Therefore, one may conclude that the wire-temperature 
traces do not s\jggest any sort of two-stage shock wave, but neither do 
they preclude the possibility of its existence. 

Comparison of wire-temperature data with optical -reflectivity data .- 
As may be seen in figure I 3 , there was unfortunately little Mach number 
overlap between the wire -temperature data and the data of Greene, Cowan, 
and Homig, and where the overlap occurred, at about M^_ = 2, the two 

types of experimental shock thickness are in rather poor agreement. How- 
ever, of the optical-reflectivity data, only the relatively isolated 
points at =2.1 differ very significantly from the Navler-Stokes pre- 
dictions. The rest of these data may be fitted quite satisfactorily by 
Navier-Stokes theories, with bulk viscosity equal to zero for the mona- 
tomic gas and with a bulk viscosity, for the diatomic gases, approxi- 
mately eqiial to the value used to fit the wire -temperature data for air 
and nitrogen. Because of the extreme difference in the experimental 
methods, this implied agreement between the present results and those 
obtained by Greene, Cowan, and Homig seems quite remarkable. 


CONCLUDING REMAEKS 


The profiles and thicknesses of normal shock waves of moderate 
strength have been determined experimentally in terms of the variation 
of the equilibrium temperatirre of an insixlated transverse cylinder in 
free -mole exile flow. The shock waves were produced in a steady state in 
the jet of a low-density wind txinnel, at initial Mach nximbers of I.72 
and 1.82 in helixjm and I.78, I.85, 1.90, 1.98, 3*70, and 3.91 in air. 
The shock thickness, determined from the maximxim slope of the cylinder 


temperature profile, varied from 5 to 3^ times the length of the Maxwell 


mean free path in the supersonic stream. A comparison between the exper- 
imental shock profiles and various theoretical predictions leads to the 
tentative conclusions that; (l) The Navier-Stokes equations are adequate 
for the description of the shock transition for initial Mach mmibers up 
to 2, and (2) the effects of rotational relaxation times in air can be 
accoxmted for by the introduction of a "second" or "bul k" viscosity coef- 
ficient equal to about two-thirds of the ordinary shear viscosity. 


University of California, 

Berkeley, Calif., May 2k, 195^ • 
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APPENDIX A 


NUMERICAL EVALUATION OF CONDUCTION AND RADIATION LOSSES 
FROM EQUILIBRIUM-TEMPERATURE CYLINDER 


A very comprelieiisive study of tlie effects of conduction and radi- 
ation heat losses on the temperature of a hare thermocouple or resistance 
thermometer in contact with a gas has been given in reference 38. This 
has been used, in conjunction with empirical information supplied by 
the authors of reference 29, to estimate the severity of these effects 
in the present experiments. 

To begin with, the observation made by the authors of reference 29 
in experiments involving the free-molecule cylinder in a \miform stream 
was that radiation losses from all wires tested were practically negli- 
gible but that end losses from their models were serious. The first 
conclusion agrees nicely with information on the radiant emissivities 
of iron and tungsten (ref. 39), giving them to be all about 0.1. This 
may be seen by evaluating the relative magnitudes of the radiation and 
convection terms in the full heat-balance equation (ref. 29): 


(j + Jj-)f(s) ^ = a(s) 


pVma 


€d(Tw^ 



(Al) 


As the worst conditions encountered in the present tests, take 
j = 2 (in a(s), as in eqs. (B24)) 
s = 1.4 

p = 38 n Hg = 0.106 Ib/sq ft (static pressure) 

Vjn = 1,074 ft /sec (most probable molec\lLar speed) 

a = 0.9 (accommodation coefficient) (refs. 29 and 40) 
e = 0.1 (radiant emissivity) 

T^ = 598° R 

Ts = 523° R (temperature of radiant surroundings) 

a = 3.74 X 10“^*^ ft-lb/(sq ft )(°R)^(sec) (Stefan- Boltzmann constant) 
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The first term on the right of equation (A1) represents convection and 
has the magnitude 60.2. The second or radiation term has the magni- 
tude thus causing a O.7 percent or about 4.3° F drop in T^. 

It is interesting to make a companion calculation of the radiation 
loss downstream of a normal shock wave at this same flow setting. There, 

s = 0.536 

p = 0.340 Ib/sq. ft 
= 1,296 ft /sec 


Tw = 554° R 


(All other quantities are the same as in the above example.) In this 
case the convective term has the magnitude 2 k . and the radiant term 
amounts to 0.0^0, causing a O.I7 percent or 0.9^ F drop in Tv It is 
seen that in this worst case radiation losses can have a measurable but 
unimportant effect on the shock profile. 

In order to evaluate the effects of end conduction losses on the 
temperature distribution along the cylinder^ an approximate value of the 
heat-transfer coefficient due to free-molecule convection is needed first. 
This is given in reference 29 in a form equivalent to 


he - 


j + ^ 


f (s) 


T 


(A2) 


This reaches its minimum value, within the present experimental range, 
at the flow setting cited above. Ahead of the shock, where T = 356° R, 
its value is 


^ ^ (0.106) (1,074) (0.9) (16.96) 

(ii. 12)(336) 

= 0 ,k 66 ft-lb/(sq ft)(sec)(*^) 


= 2.2 Btu/(sq ft)(hr)(°R) 
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Downstream of the shock, where T = 490° R, 


(o.34o)(1.296)(o.9)(i0.72) 3.600 


(11.12) (490) 778 


= 3.6 Btu/(sq ft)(hr)(°R) 


To continue with the end-loss calculation, two formulas are now 
drawn from reference 38. Both deal with the temperature distribution 
in a cylinder of length 2L, whose ends are in perfect thermal contact 
with infinite heat reservoirs at temperature Te and which exchanges 
heat with a surrounding gas, the heat-transfer coefficient being con- 
stant along the wire length. Introduce the symbols 

h heat-transfer coefficient to gas 

K thermal conductivity of cylinder 

d cylinder diameter 

Tv equilibriimi temperature of cylinder 


The first formiala gives the center-point temperat\ire of a butt-welded 
thermocouple made of wires of the same diameter but of different con- 
ductivities. Let the subscripts 1 and 2 distingxiish between prop- 
erties of the two wires, but assume h to be the same for both. Then 


T ' 
■‘•w 


equilibria temperatirre of cylinder in absence of 
end conduction losses 


and 



Ty - Ty' 
” ^w ' 



(A3) 
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This relation is plotted versus ni2_L in figure I8 for the interesting 

case where material 1 is iron and material 2 is constantan. For this 
purpose the conductivity of iron was taken to be exactly three times 
that of constantan. 

The second formula gives the average temperature of the centrally 
located fraction 2bL of the total length of a monometallic wire, 

T - T„t tanh (mL) cosh fm(l - b)Ll - sinh |m(l - b)Li 

!= =L (a 4) 

Te - T^. bmL 

which is plotted versus mL in figure I8 for the case b = 0.5* Shown 
for comparison is the case for b = 0, which simply gives the center- 
point temperature of the wire. 

A pair of conservative numerical examples will now be worked out, 
for consideration of the present data, using dimensions and properties 
of the actixal thermocouples and resistance wires. Clearly, the end-loss 
effect w in be greatest when the convective -heat-transfer coefficient 
is smallest, so the previously calc\ilated minimimi value h = 2.2 is 
used. Proceeding first to the case of an iron-constantan thermocouple 
0.002 inch in diameter, the following values are used: 

K(iron) = 55 Btu/(hr)(ft) (°R) 

h = 2.2 Btu/(hr)(sq ft)(°R) 

d = 1.67 X 10“^ ft 


So 


m(iron) = \ / ^ ~ 38.8 

(35)(1.67) 

The choice of values for L and for Tg is less straightforward, for 

several reasons. The actual physical length of the wires was about 
5 inches, but the length along which the flow properties and heat-transfer 
coefficient were constant was only about I.5 inches. Whereas the tem- 
perature at the ends of the 5-inch length is probably very close to the 
temperature of the wind-tunnel shell, it is very hard to guess an effec- 
tive end temperature to apply at the ends of the length along which the 
uniform properties prevail. Traverses made with a short thermocouple 
probe such as is shown in figure 5 have demonstrated that the eqToilib- 
rium temperature first rises as the probe moves from the free stream 
into the supersonic edge of the boundary layer and then drops rapidly 
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to a value close to the recovery temperature of the nozzle wall as the 
remainder of the boundary layer is crossed. Since the only pvirpose in 
the present analysis is to see whether end conduction losses are impor- 
tant, an arbitrary and presumably pessimistic choice will be made as 
follows : 

L = 3 A in- = 0.0625 ft 
Te - Tw' = -50° R 
ml = 2.h2 

and from figure I8 


Te - Tv- 

so that 

Tv - Tv- = 3.6° R 

The second calculation will be done for the same flow conditions 
and the same guess as to the effective length and end temperature but 

now for a re si stance -wire thermometer with b = 0 . 3 ^ d = 2.08 x 10 ^ foot, 
K = 87 Btu/(hr)(ft)(°R). With these valioes, one may calculate 


mL = 0.0625 


A)(2.2) 

(87) (0.0000208) 




From figure I8 there is obtained 


^w ~ '^w' 
Te - Tv’ 


0.05 


so that the end effect is just about one quarter of that noted for the 
thermocouple. To be conservative, this is doubled to allow for the extra 
losses down the potential leads. 


In the case of the resistance wire, which necessarily carries a 
current, there is also a possible extraneous effect 6n the equi- 
librium temperature due to the electrical heating. This may be evalmted 
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by use of equation (B23) of reference 28, which gives the corresponding 
term in the energy equation to be 


PVma 


Q 


where 


Q = 


i% 

jtdL 


R being the resistance of the wire length L. For the present case, 

i = 2 X 10“'*^ amp 
R = 27 ohms 
d = 2.08 X 10"5 ft 
L = O.Oi^-2 ft 


yielding 


Q = 0.59 watt/sq ft = 0.29 ft-lb/(sec)(sq ft) 

By comparison with the first calculations of this section, concerning 
radiation losses, it is seen that this effect opposes that of radiation 
and amounts to about one- seventh and two-fifths of the latter before 
and behind the shock, respectively. 
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APPENDIX B 

THEORETICAL CALCULATIONS OF WIRE TEMPERATURE PROFILES 


Since the only known method of calculating the equilibrium tempera- 
ture of the free-molecule cylinder involves direct use of the molecular 
distribution function, the first step in the prediction of the wire tem- 
perature profile of a shock is the calculation of the distribution func- 
tion at each point within the shock. The methods by which this step is 
taken are either selected from general attacks on the theory of nonuniform 
gases, such as the Enskog-Chapman method or the thirteen-moment method, 
or specially formiilated to describe the shock-wave transition, as in the 
method of Mott-Smith. The starting point is unanimously taken to be the 
Boltzmann integrodifferential equation, which eq-uates the rates of change 
in the distribution function due to the processes of drift and collision. 
Each separate method assumes the distribution function in a more or less 
distinct form, in which some of the observable properties of the mean 
gas motion appear as parameters and which is capable of assuming the form 
for Maxwellian eqiiilibrium with streaming velocity when the observable 
properties take on their values appropriate to the uniform flows upstream 
and downstream of the shock. This assumed distribution function is then 
substituted into the Maxwell transport eqioations for the summational 
invariants of an intermole c\lLar collision (mass, momentum, and total 
molecular energy), and possibly for various other quantities as needed, 
to provide a set of equations which can be solved subject to the shock- 
wave boundary conditions to yield the spatial distribution of the observ- 
able parameters of the distribution fimction. (This description is his- 
torically incorrect in the case of the Enskog-Chapman method but is ade- 
quate for the present task.) These equations in mean quantities are 
simply the Navier -Stokes, Burnett, or thirteen-moment equations which 
have served as the starting point for many theoretical papers on shock- 
wave structiore or, in the Mott-Smith calculations, they amoiant simply 
to the Rankine-Hugoniot equations relating the shock-wave end conditions, 
plus one specialized transport equation. 


Enskog-Chapman Method, Navier-Stokes Approximation 

For the calculation using the Enskog-Chapman method, the distri- 
bution f\mction is taken in a form capable of describing both monatomic 
gases and a limited class of diatomic gases for which the exchange of 
energy between the translational and internal degrees of freedom of the 
molecule is easy. This form of the distribution function was suggested 
by reference 57 t 
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Sr 


a 


2 8uAmC^C3 

5 2nk%:^ 


_!L_ 

rikT Sr^y2kT 



^tr ^ 


Sr 


a 


mC 


a 


rik;%^ 


/mC^ . 

■ 1^, • 

.int ^ 

“Ca 1 

f2E^ \ 

\5kT ' 

\ 

• X 1 - 

/ 

5;; 

ro 

1 

1 

a 


(Bl) 


The mean flow properties which appear as parameters in this expression 
are n, u, and T. The number density of molecules 'in coordinate 
space n and the mean flow velocity u are defined in the us\ml way 
by the following averages over the distribution function; 


fi(bbEi) 

i -«> 

c» 

i 

The temperature T is proportional to the mean total random molecular 
energy and is given by the expression 

" “ * Ki)n(trr,Eb 

The mean flow velocity enters the distribution function both explicitly, 
in the terms of the form and implicitly, in the quantity C^, 


since 
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C = I - u 

C = (^y “ ■*■ (^z " '^z) 


(B2) 


where | is the total absolute velocity of a molecule, C is the random 
"thermal agitation" velocity, and |^, ^ , and the Cartesian 

components of the vector 5 . 

The distribution function, as written above, purports to give the 
number density, in a 6 + j dimensional phase space, of molecules having 

position between r and "r + dr, velocity between | and | + d|, and 
internal energy between Ei and Ei + dEi, the internal energy being 
the result of molecular rotations or vibrations. Since the present exper- 
iments are concerned primarily with well-excited molecular rotations as 
a soiorce of internal energy, attention has been restricted to a molecular 
model whose internal energy may be ass\jmed to be of quadratic form in 
some continuous variable having the range -» to » in j degrees of free- 
dom. It has been assumed, furthermore, that the translational and inter- 
nal motions of the molecule are independent, so that integrations over 
translational velocities and over internal energies may be freely inter- 
changed in order. The assumption of the form of the internal energy 
allows the summations over internal-energy states to be replaced by inte- 
grations and leads to the following useful formulas: 



(B5) 
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where 

i kT 

The distribution function (B1) contains other symbols and some 
notations vhich have not been used previously in this report and which 
are consequently listed below: 


k 

Boltzmann constant 

m 

mass of a molecule 


coefficient of shear viscosity 

K 

coefficient of hulk viscosity 


heat-conduction coefficient for random translational energy 

^int 

heat-conduction coefficient for random internal energy 


tensor indices, used with convention of simimation over 
repeated indices 

5a3 

Kronecker delta function; 1 when a = p and 0 when a ^ 3 


The first step in the theoretical prediction of the wire temperature 
profile of a shock wave, the determination of the spatial variation of 
the parameters n, u, and T, is taken by substituting expression (b 1) 
into the Maxwell transport equations (refs. I 7 and I 8 ) for the quan- 
tities in, and + Ej_, where . All integra- 

tions and summations may be easily performed, and the procedure yields the 
familiar conservation equations known as the Navier -Stokes equations for 
a compressible, viscous, heat-conducting gas. These equations are spe- 
cialized here for the case of a one -dimensional, steady flow: 
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(nmu) = 0 
dx 


dx 


nmu^ + rddT H g 


= 0 


(Blf) 


/nmul i-^ — + I I + u 
dxl \ 2 m 2 


- (5 1" “ji: 


, dT 
- ^ dJ 


= 0 


in which A = + A'tr = Total heat conductivity. 


Each of these equations may he integrated once immediately, yielding, 
with some regrouping of terms in the third equation. 


nmu = M 

nmu^ + nkT - p + ^ 

I i«u5 + I (J + 5)unKT . ^ ^ - X f = 


(B5) 


Here M, P, and Q are constants of the flow (note that this M is not 
the Mach number), which will he used in forming the following list of 
dimensionless variables (ref. I5) ; 


V = Mu/P 
T = M^kl/p^m 
r = nmPy^M^ 
a = Mq/p^ 


(b 6 ) 


the two reference lengths 
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In terms of these variables, equations (B5) may be written as 
rv = 1 


rv^ + rr - 

3 dx 


(B8) 


rv5 + (j + 5)rTv - — Lv — - i^L' — = a 

3 dx 2 dx 


The first of these may be used to eliminate r in the other two, idiich, 
when reairranged slightly, may be written as 


itL^ = v + l- l 

3 dx V 


^ L' = 2v - v^ + (j + 3 )t - a 

2 dx 


(B9) 


The boundary conditions subject to which eqioationB (B9) must be solved 
for a shock wave are 


(BIO) 

Imposition of conditions (BIO) upon equations (B9) yields simple simiil- 
taneous algebraic equations which may be solved to display the Rankine- 
Hugoniot equations in terms of the present variables. The resialts are 


^ = ^ = 0 at x=±co 
dx dx 


v = i= 1 


2(J + k) 


(j + 5) ± G 


T = 


Hi + 


j2 + 8j + 15 ^ 2e - 


(bu) 


in which the upper sign pertains to the upstream limit and the lower 
sign, to the downstream limit, and in which there appears a new shock- 
strength parameter 


56 


NACA TN 3298 


€ 


= \lii + 5)^- 4(J + *t)<i 


(B 12 ) 


Which is related to the initial Mach nvunber by 



(BI 3 ) 


€ 


(j^ + 8j + 15 )(m3_^ - 1 ) 


(j + 3) + (J + 5)Mi' 


2 


The methods of reference 15 will be continued here with the intro- 
duction of normalized velocity and temperature variables w and t, 
which run from 1 to -1 and from -1 to 1, respectively, as the shock is 
traversed from upstream to downstream. These new variables are defined 
by the equations 


If equations (B9) are rewritten in terms of the normalized variables, 
and if the second equation is then divided by the first, the distance 
variable x is eliminated and the direction field equation for the 
velocity-temperature integral is obtained in the form 





(Bl4) 




(BI 5 ) 


dw 


2(t + w) - e(l - w^) 


in which the dimensionless parameter 
number by 


is proportional to the Prandtl 
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or 



(BI6) 


For the purpose of the present calculations, equation (BI5) niay be 
integrated numerically, with the saddle-point singiilarity (refs. 15 
and 16) w = -1, t = 1 as initial poiqt and with initial slope deter- 
mined from equation (BI5) by L' Hospital's riole. For the particular inte- 
grations which were carried out, it appeared entirely satisfactory to 
taJce increments of Aw = 0.1, with an occasional, need for finer steps 
as the upstream limit was approached. The Prandtl niomber was usually 
considered constant throughout the shock, although no essential difficialty 
and little additional work is involved in using the empirical variation 
of Prandtl n^mlbe^ with temperatiore for the particular gas concerned. 

Having found t as a function of w within the shock, w may 
now be calculated as a function of x by a second numerical integration, 
this time of the moment\jm equation (first of eqs. (B9)), which now takes 
the form 


Before this is done, it is interesting to introduce a dimensionless dis- 
tance variable y, defined by 


it L — = + v) - e(l - w^) 

5 dx (j + 5) + ew 


(BI7) 



(BI8) 
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in which the reference length L* is to be evaluated at a temperature T* 
which has here been chosen as 

T* = To (BI9) 

Tq being the reservoir temperature associated with the uniform flow at 

either end of the shock. This choice is entirely arbitrairy and is con- 
sidered appropriate to the present investigation because this value of 
the temperature is reached at or near the point of maxim\jm shear stress 
and because it is very nearly a constant for each gas tested in the wind 
tunnel, which operates at essentially constant Tq. 

The choice of y as a dimensionless distance is such that the 
limiting profile for a shock of zero strength has a particularly simple 
form in terms of w and y. The form is independent of j, Pr, the 
temperature dependence of viscosity, the particular dependent variable 
(i.e., velocity, temperature, or density), and the scheme of calculation 
(i.e., Navier-Stokes or thirteen-moment) and is equivalent to Taylor's 
early solution for the profile of a vanishingly weak shock. Thus, as e 
tends toward zero, the shock profile becomes 


= I rtf 


(B20) 


As a final step previous to the next integration, equation (BI 7 ) 
may be rewritten in the form 


dy = 


+ if ^ 




(j + 5) + ew 


J + 5, 


'2 + (j 3)^i . I^Pr + “) - - '^) 


dw 


(B21) 


The ratio may be calculated from tabulated viscosity data 

for the gas and temperatures involved, or it may be approximated by some 
simple empirical fitting formula such as the power law 


i 


(B22) 
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The choice of zero point for y is entirely arbitrary; in the present 
calculations y was usually taken to be zero in the vicinity of the 
maximum dw/dy, and the numerical integration proceeded upstream and 
downstream from that point by use of the trapezoidal rule, which seemed 
sufficiently accurate for this application. 

The spatial variation of the parameters of the distribution function 
having thus been found, the calculation of the equilibriimi temperature of 
the free-molecvile cylinder proceeds by direct application and extension 
of the analysis of reference 30 * The distribution function used here is 
somewhat different from that employed in reference ^ 0 , but the arguments 
and analysis are identical and will not be repeated here. The final 
equation for the steady-state energy balance on an insulated cylinder, 
with equal accommodation of translational and internal molecular energies 
to the temperature of the wire, can be put in the form 


„ a(8) + b(6)(i!: 

(j + 4) ^ 

T 


+ c(s) 


nkT dx 


^J + d(s) 


niikT dx 




^,nukT dx ^ 


(B25) 


du\ ,/ K du\ . dT\ ^ 

f (s) + g(sV + h(B) + l(s) +0 

' 5 nkT dx/ InkT dxi InukT dx/ 


The new variables which appear are T^, the wire temperature, and s, 
the molecular speed ratio. The coefficient functions are: 



(B24) 


ko 
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where 



are the modified Bessel fionctions of the first kind, of zeroth and first 
orders. The speed ratio s is known at each point in the shock, being 
defined by 


s 



(B25) 


Equation (B25) suffices for the calculation of the wire temperature pro- 
file of the shock, since all quantities except Tw are known. For pur- 
poses of comparison with the velocity profile, or with experimental wire 
temperature profiles. Tv is made dimensionless and is normalized to 
run from 1 to -1 by the definition 


tv - 


T 


w 


T- 


w^ 


T 




T. 


- 1 


Wc 


(B26) 


The subscripts 1 and 2 refer to upstream and downstream limit values 
here, as thro\ighout the rest of this report. 

This completes the prediction of the wire temperature profile of a 
normal shock wave, according to the Navier-Stokes equations. Curves 
of w and of tv versus y have been calculated for a number of con- 
ditions, as noted in table II. 


It may be interesting to note some numerical magnitudes of the con- 
tributions to the wire temperature due to the non-Maxwellian terms of 
the distribution function. For this purpose, use may be made of the 
approximate formiila 


where T^° is the wire temperatirre in a gas with Maxwellian distribu- 
tion function and K^_ and K 2 are complicated functions of the speed 
ratio s. The second and third terms on the right are the fractional con- 
tributions to the wire temperature due, respectively, to stress and heat 


(s) 




nkTu 


^ 1 + kJs) ^ + Kj 


w 


nkT 
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flux. These terms reach a maximum absolute value in the vicinity of the 
center of the shock, at which point a few typical values are as follows: 


For helium at M]_ = l,6l (Navier-Stokes) : 


Ki(s) ^ 

’ nkT 


= 0.021 


max 


K2 (s) 


nkTu 


= -0.056 


max 


For helium at = 1.82 (Navier-Stokes): 




= 0.052 


max 


KgCs) 


nkTu 


= -0.05l^■ 


max 


For air at = I.89 fNavier -Stokes, ~ T 


K.(s) 


nkT 


5 

= 0.0l<-l 


-max 


KgCs) 


nkTu 


= -0.056 


max 


For air at M]_ = 5»68 


^Navier-Stokes, k = j n 


Kn(s) 


Pxx 

nkT 


= 0.159 


max 


Kp(s) 


nkTu 


= -0.157 


Jmax 
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Evidently the effects of stress and heat flux are mutually 

opposing, so that the combined effect upon the wire temperature frequently 
linderstates the significance of the individual terms. 

It is interesting to note the necessity for guessing a value 
of for air in order to calculate the tv profile, whereas no 

such need arose in the calculation of the velocity profile. This occ\ars 
because the total heat flux enters the energy equation for the gas alone, 
while the fluxes of random translational and internal moleculax energies 
play quite separate roles in the energy exchange between the gas and 
the cylinder. The guess made is entirely an uninformed one and amounts 

to taking the heat-conduction coefficients A^^ and A^^^ in the same 
proportion as the average random energy contents 5kT/2 and jkT. It 
is fortunate that the calculated value of tv is not very sensitive to 
the exact ratio of the conductivities. 

Another interesting observation drawn from these calculations is 
that, up to a Mach number of 2, the dominant terms in equation (B25) 
are a(s) and f(s), the terms which would result if the wire -temperature 
calculation were based simply upon a distribution function of Maxwellian 
form with parameters n(x) and T(x), whose spatial variation is deter- 
mined by the Navier-Stokes equations. The contributions of the terms 
involving velocity and temperature gradients in equation (B25) become 
increasingly important as the shock strength increases and are definitely 
significant within the Mach number range covered experimentally. 


Thirteen-Moment Method 

The distribution fiinction is written in the form (ref. I8) 


f = n 


2jtkT 


3/2 


m \ -mC^/21cT 


, . Pap ^ , 2a 


(B2T) 


The additional parameters which distinguish this expression from the 
corresponding one in the Enskog-Chapman formulation are Pap^ 4he diver- 
genceless part of the stress tensor, and q^, the heat-flux-vector com- 
ponent. These new variables are considered to be independent, rather 
than expressible in terms of velocity and temperature gradients, and their 
presence requires two extra eqmtions in addition to those expressing con- 
servation of mass, momentum, and energy. These extra equations are 
obtained from the Maxwell transport equations for the quantities Co^Cp 

and by the methods described in reference 18. The symbol p is 
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introduced as an abbreviation for nJcT and may be identified with the 
thermodynamic pressure and the mean normal stress at a point, since this 
entire representation is intended only for a monatomic gas. 

The prediction of the normal- shock velocity and temperature profiles 
by the thirteen-moment method has been carried through in reference 15, 
so that the analysis need not be repeated here. The only difference 
between Grad's calculations and those done for the present report lies 
in a slightly different choice for the reference temperature T*. 

The calculation of the wire equilibrium temperature follows the 
same lines as in the last section, being in this case quite identical 
to the analysis of reference 50. The steady-state energy -balance equation 
analogous to equation (B25) Is 


m a(s) + b(s) + d(s) — 
k IK = P ^ pu 

f(s) + g(s) + i(s) — 

V/ B\/p '■'pu 

The coefficient functions a(s), b(s), and so forth are the same as 

those listed in equations (B24) for the case j = 0. 

This method was used to calciilate the t^ eind w profiles of Just 
one shock wave, for helium at M]_ = I.6I, e = 1-5^ = 2/5, and 

p oc ipO.665^ ipj^e comparison with the corresponding Navier-Stokes pre- 
diction is shown in figure 16. 

Mott -Smith Method 

According to the Mott-Smith method, the gas within the shock is 
imagined to be a mixture of two families of molecioles, one having the 
mean velocity and temperature of the gas upstream of the shock and the 
other having the mean velocity and temperature of the gas downstream. 

The proportions of the mixture are allowed to vary with distance through 
the shock, in a manner satisfying the shock-wave boundary conditions 
plus a rather arbitrarily chosen transport equation. The distribution 
function is written as 



1^4 
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“I 



5/2 " 2kT2 


m 



(B29) 


e 


The method used to determine n]_(x) and n 2 (x) is amply described in 

reference 15. Substitution of distribution function (B29) into the 
Maxwell transport equations for mass, momenta, and energy yields the 
Rankine-Hugoniot equations relating U 2 and T 2 to u^^ and The 

variation of n with x is determined by solving a Maxwell transport 
equation for some quantity which is not a summational invariant in the 
collision integral. Mott-Smith chooses, for convenience, the quan- 
tity alternatively dnd evaluates the collision integral 

for rigid sphere molecules and for Sutherland molecules. Whatever the 
choice of the quantity used in the extra transport equation, the solution 
for the shock-wave velocity profile can be put in the form 


where B is a fimction of the initial Mach number and the viscosity- 
temperature dependence, as well as of the choice of transport variable. 
Mott-Smith feels that the solution obtained is sufficiently Independent 
of the latter choice to be a significant approximation to a solution of 
the complete Bolt 2 mann equation. 

The calcialation of t^ from distribution function (B29) has not 
been performed previously but is a very simple extension of the calc\ila- 
tion for the wire in a uniform stream. The steady-state energy-balance 
equation becomes 



(B50) 




(B51) 
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The coefficients a and f axe again the same functions as listed in 
equations (B24). The equation on the preceding page may be solved 
for Tw, from which the dimensionless, normalized variable tv may be 
formed : 



(B52) 


The transformation from the dimensionless abscissa Bx/z to the 
variable y used in the present paper is performed as follows for the 
Sutherland molecular model, which is considerably more appropriate to 
the present purposes than is the rigid sphere model: ^ Mott-Smith's 

definition. 


J — 2 "^l = X = -c 

From reference 17^ page the viscosity according to the Sutherland 

model is 


^ ^ / mkT 

16 o 2 V « T + 


Thus, the ratio of the reference length to Mott-Smith's reference 

length Z is 



16 



J + 5; 


J + 3\ 1 
J + 


1 

J~ + 4 


T^+ S 


(B33) 


The new symbol S represents the Sutherland constant for the particular 
gas of interest and is determined by best fit of the experimental vis- 
cosity data. Comparison of equations (B30), (B3I), and (B18) yields 
the Mott-Smith shock profile in terms of w and y: 



(B3'+) 


y = 


1+6 
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The function is tabulated for the transport and the 

Sutherland model in reference 13 . The introduction of viscosity coef- 
ficients and Prandtl nxmber into the shock-profile specification in 
this case is only an artificial device allowing comparison with shock 
profiles calculated by other methods, since these macroscopic gas prop- 
erties do not enter the Mott-Smith method. 

This method has been applied here to the calculation of shocks in 
helium at = I.6I (e = 1 . 5 ) and = 1.82 (e = 1 . 77 ) and in air 
at Ml = 1.89 (e = 5 ) and Mi = 5.70 (e = l+.lfO), the results being 
shown in figures 15(a) to 15(c) and 16. 
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APPENDIX C 

EQUILIBRIUM TEMPERATURE OF A FREE-MOLECULE CYLINDER 
IN A DIATCMIC GAS WITH ROTATIONAL RELAXATION 


Since the equilibrium temperature of a free-molecule cylinder is 
considerably higher in a monatomic gas stream than in a diatomic gas with 
the same velocity and temperature, the following calculation has been 
carried out to predict the temperatiore of the cylinder in a diatomic gas 
whose inert degrees of freedom are out of equilibrium with the trans- 
lational degrees of freedom. An application of the analysis is made to 
the hypothetical case of a cylinder placed Just downstream of the first 
zone of a normal shock wave followed by a rotational relaxation zone of 
comparatively infinite length. This is the shock whose properties have 
been analyzed in references 20, 21, and 22. 

The distribution fimction for the relaxing diatomic molecules is 
taken to be Maxwellian in form, and it is assijmed that it may be factored 
into a translational term and a rotational term as follows: 


f = 


' - f 


^ g-ERs/^ 


(Cl) 


where E^g is the rotational energy of a molecule in the sth rotational 
state. This is the standard approach to the free-molecule analysis of 
uniform polyatomic gas flows (ref. 28 ), the only difference being that 
it is now proposed to assign different values to the translational tem- 
perature and to the rotational temperature Tj^. In addition, separate 

accommodation coefficients for translational and rotational energies 
and aj^ are introduced, defined by: 


= ~ 

dEti - dEtw 


*^1 ~ 

<%i - <%w 
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In each of these defining equations dE represents an average energy 
flux crossing an inf initesinial section of the interface between gas and 
solid, the first subscript indicates the type of molecular energy (trans 
lational or rotational), and the second subscript has the following 
meanings : 

i pertaining to incident molecules 

r pertaining to reemitted molecules 

w pertaining to hypothetical molecules being reemitted in 

Maxwellian equilibrium at temperature of solid wall 

Since the T's and a's are talen as constants in the problem of 
determining the equilibrium temperature of a cylinder, these modified 
concepts have no effect on the integrations over the energy spaces and 
aroimd the circumference of the wire. The desired resiilts may be drawn 
from the calculations of reference 28 and will be expressed here in the 
nomenclature of appendix B: 


In this equation the molecular speed ratio is based on the translational 
temperature alone; that is, 

s = (C5) 

The speed-ratio functions a(s) and f(s) are those given in equa- 
tions (B24), with j = 0. 

Vsirious special cases obtain from different assumptions concerning 
the relative magnitude of oj. and For instance, one might argvie 

that, if the gas molecules require an enormous number of collisions among 
themselves in order to restore equipartition of energy between trans- 
lation and rotation, the rotational energy may be negligibly accanmodated 
to the temperature of an immersed solid surface, so that approximately 
Or la+ = 0. Then the energy balance would reduce to 


Es = 1 SM (C4) 

T^ 4 f(s) 

which is simply the result for a monatomic gas. 
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^9 


If, on the other hand, one subscribes to the notion of a gas mole- 
cule being momentarily absorbed by a solid surface and then reevaporated 
with little or no memory of its past history, it may appear that a col- 
lision with the cylinder is siafficiently more serious than a simple col- 
lision between gas molec\iles so that the reluctant rotational energy may 

be accommodated as well as the translational energy.^ Setting ctp = 
brings about the special energy-balance equation 


Without attempting at once to decide which of these hypotheses seems 
the more reasonable, both may be applied to the prediction of the cylinder 
equilibrium temperature in a flow where is the translational tem- 

perature immediately downstream of a shock which was too sudden for % 

to have made any appreciable adjustment away from its upstream equilibrium 
value. The result will be compared with the temperature which the cyl- 
inder would have behind a shock initiating with the same upstream con- 
ditions but displaying no rotational heat lag. 

The ratio Tj^^T^ and s for the station Immediately behind the 
first shock front must be computed. Since Tp, the equilibrium- state 
temperature upstream of the shock, is already known, the unknowns are T^ 

and u. They may be calculated algebraically by considering one- 
dimensional flow through a plane in which the specific-heats ratio of 
the gas changes from 7/5 to 5/3^ with mass, momentum, and energy being 
conserved. A convenient technique employing published gas-dynamics 
tables is as follows: 

( 1 ) Corresponding to the initial Mach number M^, calculate a 
"monatomic" Mach number Mjj^ by 


^ = 1 a(s) + 1 ^ 


(C5) 


6 f(s) 3 




ludsen (ref. 4^1) reached this conclusion for hydrogen and 


platinum. 
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(2) Look up Tta/Tti and M2m opposite M-;^ in tables of normal- 

shock functions for a monatomic gas. Then the quantities needed for the 
problem are obtained simply from 


T^\ 


T^ 


= Itl 

Tt/p 


(07) 


since %2 = ^R1 = '^tl and 


®2 



(08) 


(5) Calculate Tv^Tt2 from equation (Ci<-) or from equation (C5). 

(k) Suppose that the upstream flow has a stagnation temperatiire Tq 
so that T.^3 _^Tq can be found opposite in adiabatic-flow tables of 

a diatomic gas. Then Tv may be referred to Tq by the equation 


^ ^ 

To Tt2 Tti To 


For comparison, T^^Tq for the case with no rotational lag may be 
read from figure 11 opposite the Mach number M2, which is given as a 
function of M^^ in diatomic-gas normal-shock tables. 

In this manner the ratio Tv/Tq versus Mq has been computed and 
graphed in figiare I7 for the three cases (a) with no relaxation zone 
behind the shock, (b) with an extended relaxation zone and with zero 
accommodation of rotational molecular energy to the surface temperature 
of the wire, and (c) with relaxation and with equal accommodation of 
rotational and translational energies. 
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TABLE I 


SUMMARY OF FLOW PROPERTIES 


Assumed stagnation temperature, 530° R 


Nozzle 

Gas 

Flow rate, 
Ih/hr 

Ml 

e 

pu, 

Ib-sec/ft^ 

^1. 

ft 

(a) 

ft 

(b) 

y/x, 

ffl 

(c) 

6 

Helium 

2.88 

1.72 

1.65 

0.000312 

0.00232 

0.00108 

174 

6 

HeliiJm 

3.84 

1.82 

1.77 

.000390 

.00189 

.000864 

233 

6 

Air 

2.6 

1.78 

2.79 

.000268 

.00254 

.00184 

154 

6 

Air 

3.9 

1.85 

2.94 

.000380 

.00182 

..00150 

230 

6 

Air 

5;2 

1.90 

3.02 

.000483 

.00146 

.00102 

501 

6 

Air 

7-7 

1.98 

3.15 

.000670 

.00106 

.000736 

454 

8 

Air 

5.2 

3.70 

4.4o 

.001280 

.000515 

.000586 

1,157 

8 

Air 

10.5 

3.91 

4.46 

.002058 

.000305 

.000240 

1.886 


®Meaii free path ahead of shock, A-i = \l— (^ 

+ 3 ) pu 


^For helium, L^^ = ^x*/pu; for air, L- = — — . 

' 2 pu 

°For helium, y = — — ; for air (Bulk viscosity = Two-thirds 
55 L* 

y = 

266 L* 


shear viscosity). 
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TABLE II 


THEORETICAL SHOCK PROFILES ACCORDING TO 
NAVIER- STOKES EQUATIONS 


Gas 

Ml 

e 

j 

Pr 

k /^ 


CD^in p oc t“^ 

^Heliim 

1.61 

1.5 

0 

2/3 

0 

0 

0.663 

^Helium 

1.82 

1.78 

0 

2/3 

0 

0 

.663 

Air 

1.89 

5 

2 

3/4 

0 

2/5 

.76 

^•Air 

1.89 

3 

2 

3/4 

2/3 

2/5 

.76 

Air 

3.70 


2 

3/4 

0 

2/5 

.76 

Air 

3.70 

4.4 

2 

3/4 

4/9 

2/5 

.76 

^•Air 

3.70 

4.4 

2 

3/4 

2/3 

2/5 

.76 

Air 

5.70 

4.4 

2 

bpj.(T) 

0 

2/5 

\(T) 


p 

^•Also treated by Mott-Smith theory, transport, 

and Sutherland molecilLes. 


^^^pirical variation with temperature (refs. 5 ^ 
and 35 ), T* = R. 
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Figure 1 .- Shock holders. 
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L- 879 i^.O 

Figiire 2 .- Air-afterglow photos showing principle of shock-holder operation. 
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Figure 5*- Distvirbance of shock by small impact tube 
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Figure Shock-wave traces taken with a small Impact tube, showing 
effect of variable choking of shock holder. 
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Figure Early style thermocouple probe and no. 4 shock holder. 


ON 


MCA TN 5298 


62 


MCA TN 5298 


SIZE OF SHOCK- 
HOLDER ENTRY-^ 



(a) Resistance-wire thermometer. 



(b) Thermocouple. 

Figure 6 .- Equilibri\jm-temperat\ire probes. 
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Figure 7 *- Experimental setup showing nozzle 
no. 5 shock holder, dial gage, and wiring 
thermometer. 


L-879ii3 

6 , wire traverse, 
for resistance 


MEAN FREE PATH, IN. 
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Figure 8.- Mean free path in free stream of no. 5 wind tunnel. 
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(a) Nozzle 6 ; 2.9 Ib/hr; helium; = I.72. 
Figure 9 *- Experimental shock profiles. 
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(b) Nozzle 6; 3«8 Ib/hrj helium; = 1.82. 


Figure 9«- Continued 
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(c) Nozzle 6 ; 2.6 Ib/hr; air; Mj_ = I.78. 


Figure 9 «- Continued. 
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(d) Nozzle 6; 5*9 Ib/tu:’; air; = 1.85* 


Figure 9»- Continued 


NA.CA TN 3298 


1.0 

.8 

.6 

.4 

.2 


tw 0 


-2 

-4 

-.6 

-.8 


.0 

-3 











0 

0 RUN 483,0.00025" RESISTANCE 




WIRE, 

, NO. 4 SHOCK HOLDER 

o.oo2"thermocouple, 

NO. 5 SHOCK HOLDER 










































-2 


0 

y 


I 


(e) Nozzle 6 j 5*2 Ib/hr; alrj M]_ = I.90. 
Figure Continued. 
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Figure 9 »- Continued 
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(g) Nozzle 8j 5*2 Ib/hr; airj Mjl = 5-70. 
Figure Continued. 
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(h) Nozzle 8; 10.5 Ib/hr; air; = 3.91. 





Figure 9*" Concluded 
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(a) Maximum- slope thickness. 



(b) Area thickness. 

Figure 10.- Definitions of shock thickness. 



Figure 11 .- Comparison of experiment and theory for cylinder temperature 

in \miform stream. 
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Figure 12.- Area shock thickness versus shock strength. 
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(a) Monatomic gas. 


(b) Diatomic gas. 


Figure 13 •- Shock thicknesses 
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Figiire Maximim- slope shock thickness in terms of upstream mean free 

path. 
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(a) Monatomic gas (helium); = 1.82. 
Figure I5.- Comparison of experiment and theory 
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Figure 15 •- Continued. 
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(c) Diatomic gas (air); = 3 -JO. 


Figure 15 .- Concluded 


NACA TN 3298 



y 

35 


Figiire l 6 .- Comparison of various theories for monatomic gas (helium). 

= 1 . 6 l. 
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Figure I7.- Cylinder temperature behind a normal shock in a diatomic gas 
exhibiting rotational relaxation (appendix C) . 


NACA - Langley Field, Va. 



00 


Figure 18.- Influence of end losses on thermocouple and resistance -wire 

probes . 
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